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ABSTRACT 

Electrostatic behavior of a colhsionless plasma in the foot region of high Mach 
number perpendicular shocks is investigated through the two-dimensional linear 
analysis and electrostatic particle-in-cell (PIC) simulation. The simulations are 
double periodic and taken as a proxy for the situation in the foot. The linear 
analysis for relatively cold unmagnetized plasmas with a reflected proton beam 
shows that obliquely propagating Buneman instability is strongly excited. We 
also found that when the electron temperature is much higher than the proton 
temperature, the most unstable mode is the highly obliquely propagating ion two- 
stream instability excited through the resonance between ion plasma oscillations 
of the background protons and of the beam protons, rather than the ion acoustic 
instability that is dominant for parallel propagation. 

To investigate nonlinear behavior of the ion two-stream instability, we have 
made PIC simulations for the shock foot region in which the initial state satisfies 
the Buneman instability condition. In the first phase, electrostatic waves grow 
two-dimensionally by the Buneman instability to heat electrons. In the second 
phase, highly oblique ion two-stream instability grows to heat mainly ions. This 
result is in contrast to previous studies based on one-dimensional simulations, for 
which ion acoustic instability further heats electrons. 

The present result implies that overheating problem of electrons for shocks in 
supernova remnants is resolved by considering ion two-stream instability propa- 
gating highly obliquely to the shock normal and that multi-dimensional analysis 
is crucial to understand the particle heating and acceleration processes in shocks. 



Subject headings: supernova remnants - shock waves - plasmas - instabililes - 
cosmic rays - acceleration of particle 
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Introduction 



The discovery of tliermal and synchrotron X-rays from young supernova remnants 
(SNRs) provides the evidence that electrons are heated up to a few keV a nd that a por- 



tion of them are accelerated to highly relativistic energy in SNR shocks (iKoyama et al. 



19951 ). Because SNR shocks are coUisionless, not only particle acceleration mechanisms but 
also electron heating mechanisms in SNR shocks are not so simple. Previous studies have 
given an important key to the formation mechanism of perpendicular coUisionless shocks. 
When the Alfven Mach number Ma is larger than the critical Mach number, about 3, a 
perpendicular shock reflects some of the incident io ns to the u pstream, where a foot re- 
gion forms on a spatial scale of the ion gyroradius (ILeroyl 119831 ). The plasma in the foot 
region consists of incident ions and electrons and reflected ions and returning ions which 
are made from refle cted ions and move t o the shock after a gyration. As for the electron 
heating machanism, iPapadopoulod (119881 ) proposed that when the Mach number is larger 
than 0.5(mp/me)^/^ ~ 20, incident electrons and reflected ions excite electrostatic waves 



by the Buneman instability (iBunemaru Il958l ) because the relative velocity between them is 
large compared with the electron thermal velocity. They also suggest that after electrons 
are heated by electrostatic waves induced by the Buneman instability, ion acoustic insta- 
bility is triggered because the electron temperature becomes much higher than the proton 
temperature. As a resu lt, electrons are strongly heated by the Buneman instability and the 
ion acoustic instability. ICargill fc Papadopoulod (119881 ) performed a one-dimensional hybrid 
simulation and demonstrated that strong electron heating actually occurs. They concluded 
that an Ma = 500 shock heats electrons by a factor of 10^ across the shock. This means 
that if the upstream electron temperature is 1 eV, the downstream electron temperature 
becomes 100 keV. This value of the down stream temperatu re is much larger than the recent 



observational one for SNRs, a few keV (IStage et al. 
open issue to be resolved for a long time. 



20061 ). This discrepancy has been an 



On the other hand, IShimada fc Hoshind (120001 ) and iHoshino fc Shimadal (12002! ) per- 
formed one- dimensional full particle-in-cell (PIC) simulations to investigate the electron ac- 
celeration at perpendicular shocks. Their simulation solves a whole region of a coUisionless 
perpendicular shock and makes reflected ions self-consistently by employing a small pro- 
ton to electron mass ratio. Their results showed that electrons are not only heated at the 
foot region but also significantly accelerated by surfing acceleration mechanism. However, 
this acceleration is valid only for the one-dimensional case because the surfing accelera- 
tion strongly depends on the structure of the electrostatic potential. In our first paper 
( lOhira fc Takaharal 120071 ) . we performed two-dimensional electrostatic PIC simulations to 
solve for the two-dimensinal structure of the electrostatic potential excited by the Buneman 
instability. We employed the real mass ratio but the simulation region is limited to the foot 
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region. Our results showed that obhque modes grow as strongly as the modes parallel to the 
beam direction, that the potential structure becomes two-dimensional and that no efficient 
surfing acceleration occurs, while electron heating occurs. Thus, the problem of electron 
acceleration has been back to the start again. In that paper, we concentrated on the stage of 
the Buneman instability and did not follow the long time scale evolution after the Buneman 
nstability has saturated. 

In this paper, we study the time evolution of electrostatic collisionless plasma instabili- 
ties in the foot region by making linear analysis and by performing two-dimensional electro- 
static PIC simulation. We perform simulat ions with a higher resolut ion, a larger simulation 
box and a longer simulation time than in lOhira fc Takaharal (120071 ). Especially, we focus 
on the evolution after electrostatic waves excited by the Buneman instability have decayed. 
Our simulation substantially improves most previous works that are one-dimensional, em- 
ploy an artifically small proton electron mass ratio or impose rather strong m agnetic field. It 
is ob v ious that a rnulti-d i mensional analysis is nec essary as discussed above (IBludman et al 
1960l : iLampe et al.lll974l : lOhira fc Takaharal 120071 ). Our motivation for employing the real 
mass ratio is as follows. For a small mass ratio, the foot region in the simulation is shorter 
than the realistic one and the time scale on which electrons stay in the foot region in the 
simulation is also shorter than the realistic one because the size of the foot region is about 
the ion gyroradius rripVdc/eB, where Vd and B are the drift velocity of refiected protons and 
the magnetic field, respectively. Because refiected ions have a large free energy, we expect 
that more energy is transported to electrons through collective instabilities with the realistic 
mass ratio in the foot region. The drift velocity is not large enough to excite electromagnetic 
waves, so that electrostatic waves are more important. 

In §2 we perform linear analysis for two-dimensional electrostatic modes. In §3 we de- 
scribe the initial setting of the PIC simulations and numerical results, followed by a discussion 
in 54. 



2. Linear analysis 

In this section, we perform linear analysis for two-dimensional electrostatic modes in 
unmagnetized plasmas with beams. In the foot region of perpendicular shocks of SNRs, 
we regard that there are several beams with a finite temperature and that their relative 
velocities are much smaller than the light speed. Therefore, the fastest growing modes are 
electrostatic modes. Then, we here concentrate on the electrostatic modes. For typical 
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interstellar medium, the ratio 



^ ~ 10-3 ( AS] ( ] 



-1/2 



is relatively small, where ilce, <^pe and Ue are electron cyclotron frequency, electron plasma fre- 
quency, and electron number density, respectively. So plasma oscillations are hardly changed 
by the magnetic field in the foot region of shocks of SNRs. When we consider spatial scale 
smaller than the gyroradius, we may neglect the effects of magnetic fields. Thus, we concen- 
trate here on unmagnetized plasmas. 

We define such that the x-direction is shock normal direction and the y-direction is the 
direction that is perpendicular to shock normal and wave vectors are on the x — y plane. For 
unmagnetized coUisionless plasmas, the electrostatic dispersion relation reads as 



sO 



k- V 



0, 



(2) 
(3) 



where the subscript s represents particle species, here electrons, ions and beam ions, tUps = 
{Aixrise^ /rrisY/'^ is the plasma frequency of the particle species s and fso is the normalized 
distribution function of the particle species s. 
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where t'th.s = (2A;Ts/ms)^/^ and va,s is the thermal velocity and drift velocity of the particle 
species s, respectively. 



To make equation (2) simpler, we use new coordinates x' and y' as 

x' — X cos 6 + ysm 9, 

y' — —X sin 9 + y cos 9, 
h 

cos 9 — 



(5) 



Then, equations (2) and (3) become 
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Finally, we substitute equation (7) into equation (6), and we obtain 



2i ;2 

-^th,s 
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wher e 2'(^s) is the plasma dispersion function and it can be numerically solved (jWatanabe 



199ll). 



We here present results of the linear analysis about two cases of plasma conditions and 
discuss on three kinds of plasma instabilities. 



2.1. Buneman instability 

We first consider the situation in which there are three beams, incident protons, incident 
electrons and reflected protons and the temperatures of all plasma beams are low, typically 
around 1 eV. We thus we neglect the contribution of returning ions in the dispersion relation, 
for simplicity. As is easily understood, the returning component plays a role in assuring 
the vanishing net current in the unperturbed state. In the dispersion relation, they play 
a symmetrical role to the reflected ions and their effects are starightforwardly understood 
when we make clear the role of reflected ions. We make analyses in the upstream rest frame 
in which only reflected ions have a drift velocity, typically t>d = fd,rof = 0.02c = 2vsh, where 
Vsh is the shock velocity. Hence, a typical velocity ratio is Vd/vth,e = 10. We assume that the 
proton reflection ratio is n^ei/^p = 0.25. 

The growth rate obtained by solving the linear dispersion relation is displayed in Figure 
1(a). In this condition, the most unstable mode is the Buneman instability. The Buneman 
instability is caused by the resonance between the electron plasma oscillation of the upstream 
electrons and proton plasma oscillation of the reflected proton beam. In Figure 1, k^ and 
ky are wavenumbers normalized by uj^^/va and the color contours show the growth rate 
normalized by cUpe, where only the growth rate of growing modes is shown. As is seen, the 
growth rate of obliquely propagating modes is as large as that of modes parallel to the beam 
direction. This feature of the Buneman instability can be well understood in the cold limit. 
We present the results of the cold limit for which all temperatures are set to zero in Figure 
1(b). In the cold limit, the distribution function (4) becomes 

fso = 5{v^-Vd^^)5{vy), (9) 
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and the dispersion relation (2) is reduced to 



ps 



0. 



(10) 



Because ky does not appear in equation (10), the growth rate of electrostatic instabilities 
in the cold limit does not depend on ky. Therefore, in the cold limit, excited waves have 
any ky and the structure of electrostatic potential to the y-direction is strongly disordered 
and loses coherence to the y dir ection. This feature is ve ry important to negate the electron 



surfing acceleration mechanism (jOhira fc Takahara 



Buneman instability in the cold limit is (iBuneman 



19581 ) 



2007 ) . The maximum growth rate of the 
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In reality, because of a finite temperature, the modes with large wavenumbers are sup- 
pressed to grow. This is seen in Figure 1(a); the growth rate for large ky decreases. The 
dispersion relation of the Buneman instability depends on Vth,e/vd and the number density 
ratio. When Vt\elvA is small, modes that have a large ky can grow as long as the wavelength 
is larger than the electron Debye length. The wavenumber corresponding to the electron 
Debye length is ZcD.e = i^pe/^^th,e = lO^Bun in the present plasma condition and the boundary 
between growing and damping region in Figure 1(a) is about ky = T/cBun ~ ^D,e, as is con- 
sistent with the present consideration. This result implies that in SNRs condition, oblique 
modes propagating to the beam direction, i.e., to the shock normal, can grow as strongly as 
the modes to the parallel direction. Because these modes are electrostatic, the direction of 
the excited electric fields is to the wave vector and the energy density of the electric field of 
the ?/-component can be larger than that of the x-component. 



2.2. Ion two-stream instability and ion acoustic instability 

After the upstream electrons are heated by the Buneman instability, the thermal velocity 
of electrons increases up to Vth,e ~ ^d while the thermal velocity of protons is roughly the same 
as the initial one, so that the situation of Te ^ Tp is realized. The other parameters are the 
same as in the low temperature case described in the previous subsection. In this condition, 
while the Buneman instability is stabilized, other types of instabilit ies can occur. We found 



that in addition to the ion acoustic instability discussed previously (jCargill &: Papadopoulos 



19881 ). the ion two-stream instability that has not been well noticed in the literature becomes 



unstable. The ion two-stream instability is caused by the resonance of the ion plasma 
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oscillation of the upstream plasma and that of reflected ions in the situation Te » Tp. Thus, 
it occurs concurrently with the ion acoustic instability which is caused by the resonance 
between ion acoustic waves of the reflected protons and electron plasma oscillation of the 
upstream plasma, where the former modes are mediated by the presence of the hot upstream 
electrons. 

The numerical results of the growth rate for ion two-stream instability and ion acoustic 
instability are shown in Figures 2(a) and (b), respectively. In Figure 2, k^, ky and the 
growth rate are normalized in the same way as in Figure 1, and the color contours show the 
growth rate of the growing modes. Note that the difference of the range of the wavcnumber 
in the x-direction kx between (a) and (b). While ion acoustic instability grows at around 
kx ~ i^pe/^d; ion two-stream instabily grows at much smaller ky. It is also noted that 
ion two-stream instability is seen for low but finite values of /c^^. In the present conditions, 
the wavenumber corresponding to the ion Debye length is A;d,p = <^pi/vth,i = lO/csun while 
that to the electron Debye length is /cBun- Ion plasma oscillation exists between these two 
wavenumbers. For wavenumbers lower than fceun, it reduces to the ion acoustic mode while 
for wavenumbers higher than A;d,p it damps by thermal motions of ions. It should be noted 
that in Figure 2(a), ion two-stream modes of parallel propagation to the beam direction 
are only weakly growing, but that highly oblique modes grow very fast and the maximum 
growth rate is larger than that of the ion acoustic instability by a factor of a few. The reason 
is explained as follows. To excite the ion two-stream instability, the resonance condition, 
kxV^ ~ uj^i must be satisfied in addition to the wavenumber condition mentioned above. The 
resonance condition requires a small kx ~ {rrie/ m^y/'^k-Q^^. For propagation parallel to the 
beam direction, this is imcompatible with the wavenumber condition and the growth rate 
is very small. In contrast, when the wave has a large fc^, both the wavenumber condition 
and the resonance condition are fulfilled simultaneously and a larger growth rate is obtained. 
The results shown in Figure 2 (a) are fully consistent with this picture of the ion two-stream 
instability. 

To understand the ion two-stream instability through the dispersion relation, we consider 
a situation where the proton temperature is zero and the electron temperature is very high 
and ky ^ kx (if the electron drift velocity is much smaller than thermal velocity, we do not 
need the final condition). Then, we can approximate as |^p|, |(^ref| 3> 1 and |^e| 1, and the 
dispersion relation becomes 



The second term represents the Debye shielding effect of hot electrons and becomes small 
for k ^ A;D,e as discussed above. It is seen that for k ^ A;D,e, the dispersion relation has the 




(12) 
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same form as equation (10) and we obtain the maximum growth rate of the ion two-stream 
instabihty when ^ref <^ JT-p, by replacing cupe with cupp in equation (11), as 



Because the upstream plasma stays in foot region by about the proton gyro-period and 
because flcp/^pp ~ 4.3 x 10~^, ion two-stream instability can grow enough in the foot region. 

As far as we are aware, this oblique unstable mode has not been considered up to now. 
We expect that this instability heats ions and that much affects subsequent electron heating 
processes. 



To perform two-dimensional simulations with real proton electron mass ratio, we confine 
our attention to the foot region through a proper modeling instead of solving the whole shock 
structure. Our simulation box is taken to be at rest in the upstream frame of reference, 
i.e., that of incident protons and electrons. We do not solve electromagnetic waves and 
concentrate on electrostatic waves. 



We define the x-direction as the shock normal pointing to the shock front, and thus the 
reflected protons move in the — a:-direction and returning protons move in the x-direction. 
The magnetic field is taken to be spatially homogeneous pointing in the z-direction and 
we solve the particle motion and electric field in the x — y plane. As the initial condition, 
we prepare upstream electrons, upstream protons, rcfiected protons and returning protons. 
Each population is uniformly distributed in the x — y plane and their momentum distribution 
is given by a Maxwellian at the same temperatures T — Tg — Tp = Ty-^i = T^et = l-75eV. 
In addition, reflected and returning protons have an extra drift velocity in the x-direction 
of Vd — ±0.04c {vd — '^Vsh)- The number densities of each population are taken as rie = 
l.Bup — lcm~^ and ^ref = n^et — 0.25np, where subscripts e, p, ref and ret represent 
upstream electrons, upstream protons, reflected protons and returning protons, respectively 
(see Figure 3). These parameters are typical of young SNRs and satisfy the charge neutrality 
and a vanishing current. 




(13) 



3. 



Simulation 



3.1. Setting 



We employ the periodic boundary condition both in the x- and y-directions. The electric 
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field is solved by the Poisson equation. We have examined two cases of the background 
magnetic field, and 90 /iG {VLcc/^pe = and 0.03), where VL^e = eB /rrieC is the electron 
cyclotron frequency. We sometimes refer the former and latter cases to the unmagnetized 
and magnetized cases, respectively. 

The size of the simulation box to the x- and y-directions is taken to be Lx = 64ABun and 
Ly = IGABuri) with a total of 2048 x 512 cells, where Abuh = Svrcijpc/t'd is the wavelength of 
the most unstable mode of the Buneman instability. Thus, the length of each cell Ax = Ay 
is 3 times the initial electron Debye length. The number of macroparticles is taken so that 
initially each cell includes 96 electrons and 96 total protons. The time step At is taken as 
5 X 10~^uj~^ and the simulation is followed until 3 x 10^uj~^ or 1.6 x 10~^r2^^ where fl~p 
corresponds to the time scale the upstream plasma stays in the foot region. 



The differences from previous simulations (jOhira Sz Takaharall2007l ) are as follows. First 



the initial temperature is lower than the previous one 7eV. This is a more realistic one 
because the typical temperature of the interstellar matter is about leV. Secondly, we add 
returning proton beam in order that the total current vanishes, although in the electrostatic 
simulation it is not so critical. Thirdly, simulation time and simulation box are larger than 
the previous values so that we can investigate the ion two-stream instability. 



3.2. Results 

Although we have performed simulations for two cases of the magnetic field strength 
(0, 90/iG), the results turn out to be almost the same. Hence, we present the results of the 
unmagnetized case and add those of the magnetized case when necessary. 

First, we discuss the time development of the electric field. The evolution of the spatially 
averaged energy density of the electric field is shown in Figure 4. The solid and dashed curves 
show the X- and ^/-components, respectively, and bold and thin curves are unmagnetized and 
magnetized cases, respectively. In the first stage for t < 250co'p|,^, the Buneman instability 
occurs and the electric field to both directions grow. After they attain peak values, they 
continue to decay till around t ~ lO^u'^. Then, after t > 10^uj~^, the y-component of the 
electric field starts to grow again while the x-component continues to decay. This feature is 
due to the ion two-stream instability as discussed below. 

It should be noted that in the first stage of the Buneman instability (0 < t < 250^;"^^), 
the ^/-component of electric field is larger than the x-component. This is different from 
our previous result. As discussed in §2, this is because the temperature is lower and 
the waves with a larger obliqueness grow faster compared with our previous simulation 
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( lOhira fc Takaharal 120071 ) . In this stage, only electron temperature has risen up to about 
me^d, but the ion temperature little changes (see Figure 5). Here, we define the temperature 
by the velocity dispersion, Tg = ms{{v — {v)Y)/2kB- 



In contrast to Cargill fc Papadopoulos ( IQSsI ). after the electrostatic waves caused by 



the Buneman instability decay at about t = 10^uj~^, only the y-component of the electric field 
grows and oscillates after the amplitude saturates. In contrast, the a;-component continues 
to decay. Of course, this feature can not be seen in one- dimensional simulations. The growth 
of the ^/-component of electric field is caused by highly oblique ion two-stream instability. 
If we consider only the parallelly propagating modes as in the one-dimesional simulations, 
ion acoustic instability has the largest growth rate. In the two-dimensional simulations we 
can take obliquely propagating modes with ky ^ into account, and the growth rate of the 
ion two-stream instability is larger than that of the ion acoustic instability as mentioned 
§2. A snapshot of the electrostatic potential at the saturation phase of the ion two-stream 
instabilities is shown in Figure 6. In Figure 6, spatial coordinates are normalized by v^/upe, 
and the color contours show the electrostatic potential normalized by mef^/2. Typical wave 
length scale to the x-direction is the order of the resonance scale 2'7iVd/(^pi and that to the 
?/-direction is comparable to 27r//cD,e as the linear analysis predicts. Before the saturation 
stage, the wavelength is smaller than 27i electron Debye length for which the linear growth 
rate is higher. 

Time development of the temperatures is shown in Figure 5. The electron temperature 
rises rapidly up to about nievl when the Buneman instability grows, while at this phase the 
proton temperature does not rise so much. At about t = 10^u~^, the proton temperature 
begins to rise by the growth of the ion two-stream instability but the electron temperature is 
kept at almost the same because this instability occurs between two ion beams and the energy 
density of the electric field is nearly 100 times smaller than the thermal energy density of 
electrons at this stage. At the end of the simulation, upstream proton temperature becomes 
about 100 times larger than the initial value and the electron to proton temperature ratio 
Te/Tp becomes about 10. For the magnetized case, Te/Tp becomes about 7. Although it 
is not explicitly shown in this paper, the proton distribution has a large anisotropy. Only 
the proton temperature of the y-direction rises up, but that of the x-direction is kept to be 
almost constant. 

Figure 7 shows the energy distribution of electrons at the end of the simulation. The 
bold and thin curves represent non-mangetized and magnetized cases, respectively. It is 
noted that no high energy tail is seen. Even when there exists a magnetic field, electrons 
are not accelerated, that is, no surfing acceleration occurs in two-dimensional simulations. 
Both curves are of a fiat top shape and can not be fitted by a Maxwellian distribution. If 
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one draws Maxwellian di stribution in Figur e 7, it is a straight line with an inchnation of 
l/k^Tf,. As discussed in lOhira fc Takaharal (120071 ). at the end of simulation, the electron 
temperature beconaes Tp ~ O.Sme^^. These results are consistent with our previous result 
Johira fc Takaharal 120071 ) . 



4. Discussion 

Now we discuss the final outcome of the ion two-stream instability. It becomes unstable 
when two conditions are satisfied; one is the resonance condition, k^v^ ~ tUpi and the other 
is that the wave length be between the ion Debye length and the electron Debye length. 
Thus, the ion two-stream instability becomes stabilized when Tp ~ T^. At the end of our 
simulation, the ^/-component of electric field has not decayed still completely. So, we expect 
Te/Tp < 10 in the final stage, and probably ions will be heated up to Ti ~ mev\ ~ 4mefgh at 
the foot region in high Mach number perpendicular shocks. Of course, because the proton 
temperature of the drift direction is still cold in the present simulation, we must check the 
isotropilazation process of ion velocity distribution by doing longtime full PIC-simulations. 

As for the electrons, the electron temperature in the foot region is also about mefg^- 
In the later stage, the growth of the two-stream instability dominates o ver the ion acoustic 



instability and little electron heating occurs. Hence, as mentioned in lOhira fc Takahara 



(120071 ). if other electron heating mechanisms do not exist, after passing the shock front, 
electrons will undergo the adiabatic heating and finally, the electron temperature in the 
downstream becomes 

Te ~ 4 X ]^m,{2v,^f = 0.41keV (^)' , (14) 

where we assume that the compression ratio is 4. This has a very important implication for 
the electron heating process of the SNR shocks. The proton temperature in the downstream 
is Tp = 3mpt>gj^/16, hence the ratio of two temperatures is 

^ ,^ 128 rric , , 

Te/Tp ~ ^ ~ 0.023. 15 

6 nip 

This value is close to the observ ed value as long as the shock velocity Vsh is larger than 



1500km/s (lAdelsberg et al.ll2008r). Name l y, we expect that the overheating problem of elec- 



trons raised by ICargill fc Papadopoulod (119881 ) can be solved by the ion two-stream insta- 
bility. 



In this paper, we prepare three ion beams as the initial condition. We have also per- 
formed simulations for other initial conditions such that there exist upstream electrons. 
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upstream and reflected protons. Electrons and reflected protons have drift velocities to sat- 
isfy the vanishing current condition. This situation corresponds to the initial phase of the 
shock reformation phenomena. The results turn out to be basically the same as those in the 
results presented in this paper. 

Our simulation does not include electromagnetic modes. In the linear stage of the 
Buneman and ion two-tream instabilities, the effects are negligible because the growth rates 
of electromagnetic modes are smaller than those of electrostatic modes for ~ 0.01c. In the 
nonlinear stage, electromagnetic modes may be important. One of the reasons is that the 
electrostatic waves with a wave vector almost perpendicular to the drift direction may make 
currents because the charge fluctuation can not be shielded by electrons in this situation, 
where the fluctuation scale is smaller than the electron Debye length. Consequently the 
current may make the magnetic fleld. It is an interesting speculation that the magnetic 
fleld might be amplifled more rapidly by the ion two-stream instability than the ion Weibel 
instability. The other reason is that anisotropic ion heating caused by the highly oblique ion 
two-stream instability excites the Weibel instability due to ion temperature anisotropy. At 
the end of our simulations, the ion temperature of x-direction is almost the same as initial 
one and the ratio of the ion temperature of y-direction to that of x-diretion Tjy/Tix is about 
100. These two features may lead to magnetic fleld ampliflcation and accompanying particle 
acceleration and heating in the shock foot region. We will make fuU-PIC simulations in 
future work to investigate these issues. 



5. Summciry 

We performed linear analysis of two-dimensional electrostatic modes and electrostatic 
two-dimensional PIC simulations with the real proton electron mass ratio to investigate the 
time evolution of electrostatic waves in the foot region of coUisionless shocks with a high 
Mach number. We consider only the foot region by properly modehng the effects of reflected 
and returning protons. Performing the linear analysis, we have shown that after electrons 
are heated by the Buneman instability, the fastest growing mode is not the ion acoustic 
instability but the highly oblique ion two-stream instability. The latter mode, which has not 
been noticed previously, is excited by the resonace between ion plasma oscillations of the two 
proton beams when the electron temperature is much higher than the ion temperature. The 
PIC simulation confirms that the excitaion of the ion two-stream instability occurs faster than 
the ion acoustic instability and that protons are heated preferentially to the perpendicular 
direction to the shock normal direction. As a result, electron heating basically stops at 
the stage of the Buneman instability and the expected electron temperature in supernova 



re mnants is fully compatibl e with the observation, avoiding the overheating problem raised 
by ICargill fc Papadopoulod (119881 ). 
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Fig. 1. — The color contour plot of the growth rate of the Buneman instabilitiy. The left 
panel (a) is for fd/fth,c = 10, Te = Tp, while the right panel (b) is for the cold limit. 
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Fig. 2. — The color contour plot of the growth rate of the ion two-stream instability (the 
left panel (a)) and that of the ion acoustic instability, for i'd/fth,e = ^,Te = lOOTp. Note the 
diffenence in the scale of the abscissa between (a) and (b). 
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Fig. 3. — The particle distribution in the initial state for the simulation. Solid and dashed 
curves represent electron and ion distribution functions of v^-, respectively. 
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Fig. 4. — The development of the energy density of electric field. The bold and thin lines rep- 
resent non-magnetized and magnetized cases, respectively. Solid and dashed curves represent 
Ex and Ey, respectively. 



-18- 




Fig. 5. — Time development of the temperatures. The bold and thin lines represent non- 
magnetized and magnetized cases, respectively. Solid and dashed curves represent electron 
and proton temperatures, respectively. Dotted curves represent electron proton temperature 
ratio. 
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Fig. 6. — The contour of the electrostatic potential at t = 1740co'pg^.The potential is normal- 
ized by mev\/2. Note that the structure is filamentary so that the electric field is almost 
perpendicular to the shock normal direction. 
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Fig. 7. — Energy distribution function of electrons at the end of simulations. Bold and thin 
curves represent non-magnetized and magnetized cases, respectively. 



